Probing non-Gaussianities in the CMB on an incomplete sky using surrogates 
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We demonstrate the feasibility to generate surrogates by Fourier-based methods for an incomplete 
data set. This is performed for the case of a CMB analysis, where astrophysical foreground emission, 
mainly present in the Galactic plane, is a major challenge. The shuffling of the Fourier phases for 
generating surrogates is now enabled by transforming the spherical harmonics into a new set of 
basis functions that are orthonormal on the cut sky. The results show that non-Gaussianities and 
hemispherical asymmetries in the CMB as identified in several former investigations, can still be 
detected even when the complete Galactic plane {\b\ < 30°) is removed. We conclude that the 
Galactic plane cannot be the dominant source for these anomalies. The results point towards a 
violation of statistical isotropy. 

PACS numbers: 95.75.Pq, 98.70.Vc, 98.80.Es 



I. INTRODUCTION 

The search for primordial non-Gaussianities in the 
Cosmic Microwave Background (CMB) is one of the most 
important yet challenging tasks in modern cosmology. 
Any convincing detection of intrinsic non-Gaussianities 
as well as their characteristics and scaling behaviour 
would directly support or reject different models of in- 
flation, and therefore affect a fundamental part of the 
standard cosmological model. 

The currently still favoured inflationary model is 
single- field slow-roll inflation pJiSj, which should result 
in (nearly) Gaussian and isotropic temperature fluctua- 
tions of the CMB. However, deviations from Gaussianity, 
preferred directions, and other kinds of asymmetries have 
been repeatedly detected [4HT7]. already questioning the 
simplest picture of inflation. It is under discussion, if 
these asymmetries are connected to foreground influences 
p!8| 119] , which appear particularly in the direction of the 
Galactic plane. 

For investigations of CMB data sets, e.g. WMAP 
data, the analysis of Fourier phases has proven to be 
a useful method [2QH23] , since all potential higher order 
correlations, which directly point to non-Gaussianities, 
are contained in the phases and the correlations among 
them. The method of surrogate maps with shuffled 
Fourier phases [TTJ |12] represents one way of analysing 
the phases. Originally, this idea stems from the field of 
time series analysis [21H27] and describes the construc- 
tion of data sets, so-called surrogates, which are sim- 
ilar to the original, except for a few modified charac- 
teristics. The validation of these characteristics in the 
original data can then be tested by comparing them to 
the set of surrogates with appropriate measures. The 
method used in jl^ |^ tests the hypothesis that the 
coefficients aim = \ci£m\e^^^'^ of the Fourier transform 
of the temperature values T{0^(j)) have independent and 



uniform distributed phases (j)£rn ^ [— ^r, tt] calculated for 
the complete sphere S. The phases of the original map 
are shuffled, which can be done within some previously 
chosen interval of interest, A£ = [^1,^2], or simply for 
the complete range Ai = [2,imax] for some given imax- 
Every realisation of this shuffling results in a new set of 
a^rn 's, which then represents (after transforming back) 
one surrogate map. Note that every surrogate still has 
by construction exactly the same power spectrum as the 
original map. If the original map contained any phase 
correlations, these are now destroyed due to the shuffling. 
Thus, any detected differences between the original and 
a set of surrogate maps reveals higher order correlations 
and therefore deviations from Gaussianity. 

One major problem in CMB analyses is the treatment 
of the Galactic plane, which strongly influences the mi- 
crowave signal. It is possible to cut out the foreground 
affected regions [28:, but this procedure itself can affect 
the subsequent analyses. When applying a sky cut, or- 
thonormality of the spherical harmonics no longer holds 
on this new incomplete sky, which leads to a coupling 
of the a^^'s, making a naive phase shuffling impossible. 
However, one can transform the spherical harmonics into 
a new set of harmonics, which forms an orthonormal basis 
on the incomplete sky [29H3T] , where phase manipulation 
can then take place again. 

The problem of incomplete data also occurs in time 
series analysis by means of surrogates. Here, gaps can 
be overcome e.g. by the use of simulated annealing for 
reproducing the autocorrelation function [32l|33]. Still, 
the quality of surrogates constructed with this method 
seems to be questionable, since it is not ensured that no 
phase correlations are induced. 

In this work, we combine the cut sky methods with 
phase shuffling, thus enabling investigations by means of 
surrogates on an incomplete sky. Our method can also 
be extended for use on incomplete data sets in general. 









■ 




0.40 . 


0.40 




0.40 


^ 0.40 



-0.40 0.40 



-0.40 .. 0.40 



FIG. 1: Example for the method of surrogates on a cut sky: 
A simulated CMB map in its original form (upper left), af- 
ter reconstruction with cut sky harmonics corresponding to 
the sky cut \b\ < 10° and the limitation imax = 20 (upper 
right), and after two realisations of a phase shuffling of these 
harmonics for the complete range of A£ — [2, 20] (lower row). 



II. METHODS 

A. Cut Sky Surrogates 

On a complete sphere 5, an ort honor mal basis is given 
by the spherical harmonics ^^^(s) with i > 0^—£ < 
m < i and s G S. Let the number of harmonics be 
limited by some given imax ^ Now, for any map 

/(^) = Yl^Cm^ (^imyimis), V 5 G 5 with aim ^ C, 

and for any new incomplete sky iS^^^, we want to know 
the corresponding af^ and Yf^* for representing the 
map on the remaining regions of the sphere: f{s) = 
Etm «?m^/m'(5). V s e ^^^S with r/^^ being orthonor- 
mal on S^'^^ and thus aj]^ being unique. For real valued 
spherical harmonics, this was performed in [29] and [30] . 
and later on extended in [31 . The methods presented 
there can be easily adapted to complex valued spherical 
harmonics as well. 

At first, we define the vectors 

Y{s) := [ro,oW,n,oW,l^i,i(^),...,l^^_/_W]^, 

a := [aO,0, ^^1,0, 0.1,1, Oi^a.,irr^aJ^ 

containing all harmonics and coefficients with m > 0, 
respectively, of the given map on the complete sphere. 
Both have a length of imax {imax + ^){imax + 2)/2. 
Analogously, we define Y'^'^'^{s) and a^^^. Our objective 
is to determine two transformation matrices Bi^ B2 G 
£^maxx^max ^ that fulfil thc followlug cquatious: 

r^^^(5) =5i Y{s) (1) 
a^^' =B2a (2) 

To identify them, we need to define the coupling matrix 

C := [ Y{s)Y%s)dn 
Jr 
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FIG. 2: The latest KQ75-mask provided by the WMAP-team 
[28] , that cuts out the highly foreground affected regions, com- 
pared to the three central latitude cuts applied in this work. 
The smallest cut is ±10° (dark blue), which removes already 
a large amount of the highly foreground affected regions of the 
Galactic plane but presumes nearly all non-affected regions. 
The second, ±20° (blue), removes more regions, while the 
largest, ±30° (light blue) cuts out almost the entire Galactic 
Plane. Only minor point sources (turquoise) remain. 

as well as analogously its counterpart C^^^ with R be- 
ing a given region on the sphere. Hereby, F* denotes 
the hermitian transposed of Y. When working with 
a pixelised sky, one uses a sum over the pixels of R 
instead of the integral. For R = S^^^ ^ an ortho nor- 
mal set of harmonics Y[^ needs to fulfill the condition 
(jcut _ J-^^^^ with li^^^ being the unit matrix of size 
^max- We can use equation ^ on C^^^ to change this 
condition to BiCBf = li^^^. It is possible to apply dif- 
ferent matrix decompositions to obtain C = AA* with 
^ ^ ([^^maxy<^max ^ Cousequeut ly, the above equation now 
reads as {BiA){BiAy = li^^^ and offers the simple so- 
lution Bi = A~^. 

For the evaluation of let us recall that the coeffi- 
cient vector a can be expressed by 

a = [ Y{s)f{s)dn 
Js 

or, respectively, 

a^^^ = / Y''''\s)f{s)dQ . 

Jscut 

Inserting ([T]) and the expression f{s) = a^Y{s) into 
the latter leads to a^^^ = BiC^a. Now we use the 
above matrix decomposition again and obtain a*^^* = 
^i(AA*)^a = A^a. Thus, it follows that B2 = A^ . 

To obtain the af^ and Y^^* with m < 0, we make use 
of the following equations, that hold for full sky and that 
we assume to be valid also on incomplete skies: 

and 
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For the sky cuts and -^-ranges used throughout our study, 
the cut sky harmonics were tested and confirmed to be 
orthogonal. 

For constructing C = AA"" as above, one can make 
use of different matrix decomposition methods. Since 
the couphng matrix C is hermitian and can be treated as 
positive definite for low ^ by construction, a Cholesky de- 
composition is applicable. This is the easiest and fastest 
way, although numerical problems only allow usage for 
lower Irnax [ST • Another possibility is the eigendecom- 
position (ED): We obtain C = VWV*, with the columns 
of V containing the eigenvectors, and W being diago- 
nal and containing the eigenvalues of C. Because of 
the properties of the coupling matrix, these values are 
real and positive, allowing therefore a simple decompo- 
sition of W by taking the square root of every element, 
W = iyV2(i^*)i/2 Thus, we obtain A = VW^''^. Since 
C is hermitian, the ED is formally similar to a singular 
value decomposition (SVD), which is also applied in this 
paper, with the eigenvalues corresponding to the singular 
values. For both the ED and the SVD we apply a house- 
holder transformation similar to [31 to make A lower tri- 
angular. For the Cholesky decomposition, this is already 
the case by definition. Thus, due to equation ([2|, it is 
ensured that the mono- and dipole contributions of the 
underlying maps - often considered as non-cosmological 
- are kept separate from the i>2 modes. 

With the help of the new cut sky harmonics ^ we 
can now generate the surrogates on a cut sky S^'^^ as well. 
Similar to above, we shuffle the phases (\f^ of the cut sky 
coefficients a^J^, which is in this work performed for the 
full cut sky range A^^^^ = \l^tmax\' We obtain new sets 
of a^g^ 's, which are transformed back to pixel space to 
form the cut sky surrogate (CSS) maps. Figure [l] shows 
as an example a simulated CMB map, its reconstruction 
with cut sky harmonics and two corresponding cut sky 
surrogates. As we did in the case of a complete sphere, 
we now search for deviations between the original data 
as well as its surrogates. However, one has to take care 
about the above mentioned properties. While the uni- 
form distribution still holds for 0^^, the single phases in 
the sets are no longer independent from each other due to 
equation In other words, the cut sky transformation 
induces phase correlations to the underlying map. To 
account for these systematic effects, we create for each 
of the input maps Nfss = 20 full sky surrogate (FSS) 
maps as explained above, with irnax = 1024 and by shuf- 
fling the phases within A£ = [2, 1024]. By comparing the 
results of the surrogate analysis for an input map and its 
FSS, we evade systematically induced phase correlations, 
since the original map as well as the FSS are subject to 
the same cut sky transformation process. Thus, we only 
search for additional signatures possibly contained in the 
phases. 

Even though the input map might feature induced 
phase correlations, it has to be pointed out that the cut 
sky transformation is reasonable. Only with a new set of 
cut sky harmonics, analysis techniques that depend on 
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FIG. 3: Sketch of the method of surrogates on a cut sky. 



the orthonormality (like the method of surrogates) are 
feasible on an incomplete sky. The transformation into 
the new regime only induces a phase coupling to the data 
set, whose consequences can be eliminated as explained 
above. The method of cut sky surrogates itself is firmly 
constructed. 

We construct the cut sky harmonics for three different 
central latitude sky cuts, that remove \b\ < 10°, 20° and 
30° of latitude in the centre of the maps. While the small- 
est cut {\b\ < 10°) removes already a large amount of the 
highly foreground affected regions but retains nearly all 
non- affected regions, the largest {\b\ < 30°) excludes al- 
most the entire Galactic plane, with only minor point 
sources remaining. This is illustrated in figure [2] We 
choose an upper bound of imax = 20 and set airn = for 
i > imax • To check for consistency with (TTJ [12] , we ap- 
plied the cut sky formalism also to the complete sphere 
with no points excluded. To compare the different matrix 
decomposition methods, all three approaches (Cholesky, 
ED, SVD) were applied. For every sky cut, the phases of 
the coefficients a^J^ were shuffled a hundred times to gen- 
erate Ness = 100 cut sky surrogates for each input map. 
The same was done for the corresponding FSS maps. 



B. Measures for comparison 

In general, the comparison of the original data and 
its surrogate maps can be accomplished with any higher 
order statistics. For our study, we chose the scaling index 
method (SIM) [TOl [12] as well as Minkowski functionals 
[H [36] as test statistics. 

The SIM is a local measure that is able to detect 
structural characteristics of a given data set by estimat- 
ing its local scaling properties. Briefly, the tempera- 
ture anisotropics T(^, 0) are transformed to variations 
in radial direction around the sphere, therefore leading 
to a point distribution p^^ i = 1, . . . , A/^^^^, in three- 
dimensional space. Then, the weighted cumulative point 
distribution p{pi^r) is calculated for every point Pi and 
a freely chosen scaling parameter r. Since we will only 



investigate the large scales in our study, we choose the 
free parameter r to be rio = 0.25 (in the notation of 
[To]), which is appropriate for these scales. Eventually, 
the scaling indices are obtained by calculating the loga- 
rithmic derivative of p{pi^r) with respect to r. 

The three Minkowski functionals measure the be- 
haviour of a given map with respect to different threshold 
values u. The fraction of the sky where the temperature 
value is larger than v is denoted as the excursion set 
i?(z/), its smooth boundary is identified by dR{v)^ and 
da and dl describe the surface element of R{i') and the 
line element of dRiv)^ respectively. Then, we can define 
the three Minkowski functionals as 

Marea{^) = da 
JR{v) 

Mperim{^) = / dl 
JdR{u) 

Meuler{j^) = din , 

JdR{v) 



with n being the geodesic curvature of dRiv). For more 
details, we refer to [4l|36]. Eventually, we sum up over 
all thresholds with the help of the appropriate cut sky 
surrogates by means of a x^-measure 

V 

for Marea(^), Mperim(^) and Meuier{^), respectively. 

For computational reasons, the resolution of the input 
maps in the corresponding HEALPix scheme [34", "SF was 
chosen to be Nside = 256 for the scaling index analysis 
and Nside = 64 for the Minkowski functionals. By test- 
ing several subsets, we assured ourselves that the results 
are only marginally affected when choosing a lower reso- 
lution. 



C. Statistical interpretation 

The results for the different maps of both the scaling 
indices and the Minkowski functionals are then evalu- 
ated in terms of rotated hemispheres: For 768 different 
angles we rotate the underlying maps and calculate the 
cr-normalised deviations 

Sl{Y) = (r^^^ - (Y^^^)) /(JyCSS 

of the pixels included in the new upper hemisphere 
between the input map and its cut sky surrogates, 
by means of the measure Y. In our case, Y = 

{oL),(Jc^,xlrea^xlerim^xluier^ with {o) and (7^, being the 
mean and the standard deviation of the scaling index re- 
sponse a(pi,r), respectively. The result is then shown as 
colour-coded pixel, whose centre is pierced by the z-axis 
of the respective rotated reference frame (see [10]-[12]). 
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FIG. 4: The foreground residuals resulting from the subr ac- 
tion of the WMAP ILC map from the full seven-year fore- 
ground reduced coadded VW-band in its original form (left) 
and added to a Gaussian simulation (right). The latter rep- 
resents the map that is used in section pIIA| 

To separate traces of possibly intrinsic phase correlations 
from those induced by the transition to incomplete sky, 
and thus to account for systematic effects, we calculate 
the statistics 

S^{Y) = (sr^iY) - {S[''{Y))) /as.ss^y^ 

for comparing the results of Si{Y) for the original and the 
full sky surrogate maps. Figure [S] illustrates the concept 
of this statistical interpretation of the method of cut sky 
surrogates. 

An important advantage of the rotated hemispheres 
5*2 (F) is the possibility to identify not only a global dif- 
ference between the original map and its surrogates, but 
also the directions where these differences occur. Thus, 
one can make a statement about possible deviations from 
statistical isotropy. One example for such anisotropics 
are asymmetries, that have been detected in numerous 
previous analyses as mentioned above. A difference be- 
tween a S2{Y) map that detected asymmetries and a sta- 
tistical isotropic map is reflected in a difference in the 
empirical distributions: For the latter, one would ex- 
pect a histogram with a narrow range of values for 5*2 (^) 
that peaks around zero. But for the first, one expects a 
broader spectrum and possibly more than one peak, since 
on the map there are at least two areas with different val- 
ues. One measure that exactly accounts for the broad- 
ness of histograms is the discrete entropy H{S2{Y)). We 
define it as 

H{x) = - ^ Pilnpi , 

i=l 

with Pi being the empirical probability (that is the rela- 
tive frequency) of the realisations of x to occur in the bin 
z, with Ni)ins describing the number of bins in a chosen 
underlying interval. 

The entropy H{S2{Y)) reaches its maximal value for a 
uniform distribution, which would correspond to a 5*2 (F) 
map that consists of all different values of the chosen 
range, which again would imply that there is no statis- 
tical isotropy at all. The minimum of H{S2{Y)) is zero. 
This occurs for the special case of a completely isotropic 
map with a distribution where only for one bin i the em- 
pirical probability pi differs from zero. Thus, a statistical 
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FIG. 5: The a- normalised deviations ^(cra) comparing a sim- 
ulated map and its 20 full-sky surrogates for the complete 
sphere (upper left) and the three different central latitude 
sky cuts \b\ < 10° (upper right), \b\ < 20° (lower left), and 
\b\ < 30° (lower right), that were constructed by means of the 
singular value decomposition. 



anisotropic map, e.g. one with a strong asymmetry, re- 
sults in a higher value for H{S2{Y)) than an isotropic 
one. 



III. VALIDATION 

A. Underlying Data Sets 

To test the new approach, we generate two types of 
Gaussian simulations of the coadded VW-band of the 
WMAP satellite via a noise-weighted sum. The differ- 
ence between these two is simply that for the one type we 
add Gaussian noise with a particular variance, given by 
Nobs 5 for every pixel of the sphere and for the other we do 
not. The procedure is the same as in JW ^ but note that 
we now apply the more recent WMAP 7- year parameters. 
We include the simulation without the noise, since only 
there one can be absolutely sure that the results are not 
distorted by any influences coming from the input map: 
The map is created by only using a given power spectrum 
and completely random phases, and should therefore be 
resistant against phase shuffling in a statistical sense. In 
addition to these two clean simulations, we applied the 
cut sky surrogate approach to another simulated Gaus- 
sian map, to which we added typical foreground residu- 
als that are still present after the template cleaning of 
the WMAP data. With the help of this second sim- 
ulation, one can examine the impact of possible alias- 
ing effects due to the chosen imax^ that could be caused 
by strong foregrounds which are cut out. The required 
residuals were computed by subtracting the WMAP ILC 
map from the full seven-year foreground reduced coadded 
VW-band. However, when performing this subtraction, 
one has to bear two things in mind: First, not only the 
coadded VW-band but also the ILC map contains in- 
fluences from the Galactic plane, although with a lower 
amplitude. Second, the resulting map is slightly nois- 




FIG. 6: Same as figure [s] but for ^(Xarea)- Note the different 
colour scaling. 



ier compared to the original foreground cleaned WMAP 
band-wise maps. Still, both remarks should not have a 
significant impact on the intention of the test. The resid- 
uals resulting from the subraction as well as the corre- 
sponding simulation with the typical foreground pattern 
are shown in figure |4j 



B. Results 



For the clean simulated Gaussian maps, the signifi- 
cances 82(0- a) and S2{Xarea) Calculated for the rotated 
hemispheres are illustrated in the figures [5] and [6| The 
differences between the original and the FSS maps are in- 
significant (nearly always 5*2 (^) < 2.5) for the complete 
sphere and all sky cuts. This holds for the results of the 
scaling indices {Y = (a), a^) as well as for the Minkowski 
functionals (F = xlrea^xlerim^xLier)' There are only 
marginal differences when one adds the Gaussian noise to 
the map, which shows that these noise fluctuations alone 
are not responsible for non-Gaussianities or anisotropics 
on the large scales. The results of the assembled map are 
very similar as well, except for the full sky, where phase 
correlations are obviously present. These results clearly 
demonstrate the validity of the approach. Although 
phase correlations were induced by the cut sky transfor- 
mation, these are subtracted away by using S2{Y). The 
results of the assembled map show that the impact of 
aliasing effects, even when strong foregrounds are present 
in the Galactic plane, is negligible. 

Only minimal differences were detected for the three 
used matrix decomposition methods, that are likely to be 
due to the random shuffle of the phases. When going to 
larger ^-ranges or more irregular sky cuts, this technical 
part of the investigation will become more important, es- 
pecially for making the transition to the cut sky possible. 
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FIG. 7: The temperature point source mask of the WMAP 
team in its original (blue spots) and in its extended form 
(additional red regions), used to obtain the VW7p and the 
VW7pe data map, respectively. 



IV. APPLICATION TO WMAP DATA 
A. Underlying Data Sets 

For the application of the cut sky method to observa- 
tions, we make use of several data maps. Two of them 
are linear combinations of the different frequency bands 
and based on the WMAP results [37l [38]: First, the 7- 
year Internal Linear Combination map (ILC7) provided 
by the WMAP team [28 and second the 5-year needlet 
based ILC map (NILC5) [S^- In addition to these two 
data sets, we examine three maps based on the 7-year 
co-added VW-band: For the first (VW7m), the galactic 
plane and several point sources of the co-added VW-band 
map are cut out by means of the KQ75-mask [37]. These 
empty regions are then filled with Gaussian noise, whose 
mean and standard deviation correspond to the values 
of the remaining map, while preserving the spatial noise 
patterns. As above, this procedure is similar to [10 . The 
second map (VW7p) is obtained by applying the same 
method but cutting out the point sources outside the 
Galactic Plane only. The temperature point source mask 
of the WMAP team is used for this purpose [37 . Since 
some of the point sources correspond to larger sources 
that might not always be completely masked out, it can 
be advantageous to extend the previous mask with the 
help of the high-latitude regions of the WMAP process- 
ing mask [37J. This leads to the third VW-band based 
map (VW7pe) used in this paper. The regions that were 
cut out to obtain the VW7p and VW7pe maps are illus- 
trated in figure [7| For all data maps, the monopole and 
dipole were removed. 

The three maps based on the co-added VW-band map 
will be treated slightly different to the analysis process 
from above: Instead of using the usual 20 FSS maps, we 
compare the VW7m, the VW7p, and the VW7pe map 
with 20 simulations of the VW-band, where we likewise 
excluded the respective mask or point sources and filled 
it with Gaussian noise. In doing so, we avoid to detect 
deviations that are only due to the structural differences 



inside the masked regions or to the peculiar boundary 
between the masks and the rest of the map. These unin- 
tentional deviations might occur when applying the FSS 
method. However, for the VW7m map, a comparison 
with FSS was tested as well, leading to similar results as 
the approach with the masked simulations. 

The ILC7 and NILC5 data sets are maps that can also 
be used in a full sky surrogate analysis. The maps were 
nevertheless analysed in this work, on the one hand to 
check for consistency with the full sky surrogates method, 
and on the other hand to elucidate the influence of the 
different map making procedures on the phase correla- 
tions. Apparently, effects of these map making proce- 
dures should in particular be present in the Galactic 
Plane. 

For the VW7m map, a full sky surrogate analysis would 
not be applicable, since obviously the filling of the mask 
has a severe influence on the phase correlations. However, 
it could be used in a regular full sky analysis, e.g. in 
comparison with simulated CMB maps of the VW band, 
for which the same mask filling method is applied. 

The remaining two maps (VW7p and VW7pe) repre- 
sent data sets that do not attempt to remove the fore- 
ground influences of the Milky Way. These maps would 
not produce reasonable results when analysed with the 
full sky surrogates approach. Thus, they are a perfect 
example to demonstrate the usage of the method of cut 
sky surrogates. By applying the method on these maps, 
we are consistent with the statement that the removal 
of foreground affected regions is to be preferred to an 
reconstruction of these areas [40]. 



B. Results 

When evaluating the scaling index results 5'2((q^)) and 
5*2 (cTq;) for the \h\ < 10° cut for the rotated hemispheres 
of all data maps described above, we detect significant 
non-Gaussianities and an asymmetry. Both features were 
already found in corresponding full-sky analyses |TTJ [12] . 
The signal for 5*2 ((a)) becomes less clear for larger cuts. 
But for 5*2 (cTc^), a significant asymmetry with a clear 
north-south direction persists for all maps when exclud- 
ing the Galactic plane. As examples, the significances 
5*2 (ctq;) determined for the rotated hemispheres of the 
NILC5 and the VW7pe map for the different sky cuts 
are shown in the figures [S] and [9] respectively. 

The Minkowski functionals show similar results: For 
Mareai^)^ ouc dctccts similar deviations between all the 
data sets and its related full sky surrogates. The only ex- 
ception is the VW7pe map, where the asymmetries are 
less pronounced than in the other maps. For the VW7m 
map, that cuts out even larger areas of foreground af- 
fected data, a clear asymmetric behaviour is found again. 
As examples, the results for Marea{^) are illustrated for 



the NILC5 and VW7pe maps in the figures 10 and 11 



For Mperim{i^) ctud Mauler the rcsults agree for the 
full sky, but get less definite for larger cuts. However, ex- 
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FIG. 8: Same as figure [5] for tfie NILC5 map. 
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FIG. 10: Same as figure [6] for the NILC5 map. 
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FIG. 9: Same as figure [5] for the VW7pe map. 
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FIG. 11: Same as figure [6] for the VW7pe map. 



perience shows that this could hkely be due to the hmited 
amount of pixels that is present especially when apply- 
ing the \b\ < 20° or 30° cuts. The decreasing number of 
pixels on the incomplete skies leads to an increasing influ- 
ence of noise and a lower influence of the intrinsic signal. 
Due to discreetness, possible effects on the tails of the 
Minkowski functionals may vanish, which in turn might 
damp the respective x^- This especially holds for the two 
Minkowski functionals Mperim{'^) and Meuierii^) that ex- 
amine complex pixel formations and thus need enough 
data points to produce statistically reliable results. 

The detected north-south asymmetry is reflected in 
the empirical distribution functions of the result maps 
5*2 (c^o;)- We compare the results of the ILC7 and the 
NILC5 maps with 10 simulations without noise. These 
simulated maps undergo the same processes as the data 
maps. That means we constructed for each simulation 
20 FSS, and again for each FSS and the original map 100 
CSS for all different sky cuts. The results are illustrated 
in figure [T2j Clearly, the histograms of the data maps 
show a much broader spectrum than the respective out- 
put of the simulations, which underlines the detection of 
statistical anisotropics. 

This aspect is in turn quantified by the empirical en- 
tropy H{S2{cFoi))' Figure [13] presents the corresponding 
results of the twelve above mentioned histograms and ad- 
ditional ten simulations with Gaussian noise. The results 
of the data maps offer higher values than the ones of the 



simulations. Also, it is remarkable that these values are 
nearly not decreasing when moving to a larger sky cut. 
As above in section [Til B[ the differences between the two 
types of simulations are negligible. 

In summary, both the scaling indices and Minkowski 
functionals detect higher order correlations and thus non- 
Gaussianities, as well as a clear north-south asymmetry 
in all data sets even when large parts of the Galactic 
Plane are removed. Only the output of the Minkowski 
functionals of the VW7pe map shows a weaker signifi- 
cance. Nevertheless, the scaling indices detect an asym- 
metry in the VW7pe map. To investigate this particu- 
larity a bit further, we repeated the CSS analysis of the 
NILC5 map where we removed the masked regions of the 
VW7pe map. For simplicity reasons, no new FSS were 
generated. A comparison of this new NILCSpe map with 
the FSS of the NILC5 data lead to the same results as 
the analysis of the NILC5 map itself. 

The different signs in the opposite hemispheres of the 
findings for S2{cFa) that appear for all data maps reflect 
a lower (higher) variability of the scaling index results for 
the northern (southern) hemisphere of the original maps 
compared to their surrogates. A higher variability im- 
plies the presence of different types of structure. Thus, 
the lower variability in the north can also be interpreted 
as a lack of power in the northern hemisphere, which is 
in agreement with former results ^ l9l lT2| |4TJ 02] . The 
regions with large deviations from isotropy become very 
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FIG. 13: The empirical discrete entropy H{S2{cra)) corre- 
sponding to the histograms of figure |12[ The blue and red 
lines represent the ILC7 and NILC5 maps, while the results 
of the simulated Gaussian random fields with and without 
noise (10 realisations each) are illustrated by the dotted and 
dashed grey lines. 



the full sky surrogate analyses of [Hj [12], where fore- 
ground effects could in principle have played a role, and 
confirms previous findings of other analyses (e.g. |9j). In 
combination with the multitude of checks on systemat- 
ics performed for the cut sky transformation in section 



FIG. 12: The histograms of 5*2 (ctc) for the complete sphere 
and the three different central latitude sky cuts \b\ < 10°, \b\ < 
20°, and \b\ < 30° (from top to bottom), for the ILC7 (blue) 
and NILC5 (red) maps as well as for 10 simulated Gaussian 
random fields without noise (grey). 



nil and for the technique of surrogates in [TTJ [12] , the 
results point towards the conclusion that the signatures 
are of cosmological origin. This would represent a strong 
violation of the Gaussian hypothesis and of statistical 
isotropy. Both assumptions are fundamental parts of 
single-field slow-roll inflation, which is therefore rejected 
at high significance by this analysis. 

A lot of effort as well as knowledge about the cut sky 
transformation process will be necessary to be able to 
extend the method of cut sky surrogates to higher reso- 
lutions. A naive approach results in an 0(^^) scaling. As 
usual in CMB analysis, this includes both calculations 
in Fourier as well as in real space. It is important to 
remove the restrictions of a low irnax to enable investi- 
gations of higher -^-ranges, as it was already done for the 
full sky case (TTJ [12]. When going to higher ^-values, it 
is of course helpful to work with a higher resolution in 
pixel space as well. The upcoming data sets of the Planck 
satellite will offer a better resolution for both Fourier and 
pixel space. 



V. CONCLUSIONS 



broad for the larger cuts, which makes it difficult to draw 
a clear statement about a more accurate direction of the 
detected asymmetry, e.g. if it was aligned with the eclip- 
tic plane or not. 

The above results also clearly indicate that both the 
detected non-Gaussianity and asymmetry on the large 
scales cannot mainly be attributed to foreground influ- 
ences of the Galactic plane. This verifies the results of 



We demonstrated the feasibility of generating surro- 
gates by Fourier-based methods also for an incomplete 
data set. This was worked out for the case of a CMB 
analysis on an incomplete sphere. Three different con- 
stant latitude sky cuts were applied. For this purpose, 
three different cut sky transformations were calculated. 
We generated 100 cut sky surrogates for every input map, 
sky cut and matrix decomposition method, which were 
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analysed by means of scaling indices and Minkowski func- 
tionals. To remove systematic effects, a second analysis 
compared the results of the original with the ones of 20 
full sky surrogate maps for each of the input maps. For 
simulated maps, no anomalies could be detected. The 
findings for the data maps show strong signatures of non- 
Gaussianities and pronounced asymmetries, which per- 
sist even when removing larger parts of the sky. This 
confirms that the influence of the Galactic plane is not 
responsible for these deviations from Gaussianity and 



isotropy. Together with former full-sky analyses, the 
results point towards a violation of statistical isotropy. 
Similar tests with the forthcoming Planck data will yield 
more information about the origin of the detected anoma- 
lies. 

Many of the results in this paper have been ob- 
tained using HEALPix [34]. We acknowledge the use 
of LAMBDA. Support for LAMBDA is provided by the 
NASA Office of Space Science. 
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